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ABSTRACT 

To understand the chemical evolution of the Galaxy, we need to understand the contribution of 
PNe in the ^He abundance. ^He abundances have been detected only in a couple of PNe, their 
abundances are consistent with the standard models, in which ^He is produced in significant 
quantities by stars of l-2Mo. However, for all the other PNe observed to date there have been 
no detections, therefore only upper limits in their abundance can be calculated. Observations 
of the ^He+ emission line using the VLA towards the PNe IC 418, NGC 6572 and NGC 7009 
were used to obtain upper limits for their ^He abundance. Because the abundance of ^He in 
H II regions, the ISM and the pro to- solar system is much lower than what is predicted, new 
chemical models were proposed. The resulting evolution of ^He, based on stellar evolution 
models, can be consistent with the values determined in pre-solar material and the ISM if 96% 
of the population of stars with mass below ~2.5Mo has undergone enhanced ^He depletion. 
This implies that unless the combined sample of PNe that has been observed so far is very 
atypical, the solution to the "The ^He Problem" lies elsewhere. However, the results presented 
here suggest that more observations are needed in order to make a strong conclusion about 
the stellar evolution models. 
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1 INTRODUCTION 

The cosmic abundance of the ^He isotope has important implica- 
tions for many fields of astrophysics. It impacts stellar evolution, 
chemical evolution and cosmology. The present interstellar ^He 
abundance, as per all the light elements, comes from a combina- 
tion of Big Bang Nucleosynthesis (BBN) and stellar nucleosynthe- 
sis ( [Wilson & Rood|1994| . H ii regions are young objects compared 
with the age of the Galaxy, so they represent the zero-age objects. 
Their ^He abundance is the result of 13.7Gyr of Galactic chemical 
evolution. Any hydrogen-burning zone of a star which is not suf- 
ficiently hot (>7xlO^K) will produce ^He via the p-p chain. Stars 
with masses <2.5M0 are the net producers of ^He. For these stars, 
p-p burning is rapid enough to produce D in situ, and enable the 
production of ^He (D+ p ^ ^He + y). More massive stars are dom- 
inated by CNO burning; although they do produce some ^He via 
the p-p chain in their outer zones, it is not enough to offset the ^He 
destruction in the interior zones ( [Dearborn et al.|1986| ). Therefore 
PNe evolved from stars with masses <2.5Mo should be substan- 
tially enriched in ^He. PNe ^He abundances are important tests of 
stell ar evolution theory (Bania et al. 2010 ). 

. ( [T995j ) presented "The ^He Problem". According to 
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standard models of stellar nucleosynthesis there should be a ^He/H 
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abundance gradient in the Galactic Disk with the highest abun- 
dances occurring in the highly-processed inner Galaxy, the ^He/H 
abundance should grow with source metallicity and the proto solar 
^He/H value should be less than what is found in the present ISM. 

After more than 30 years of effort, observations of the abun- 
dance of^HeJn_Hr^gions_^^ with these expecta- 
tions (' Rood et al.|l979[|Bania et al.|2007] . For the ^He problem to 
be solved, it will require that the vast majority of low-mass stars 
fail to enrich the ISM. One suggestion to solve this is adding extra 
mixing in the red giant branch (RGB) stage. This extra mixing adds 
to the standard first dredge-up to modify the surface abundances. 
[Eggleton et al.| p006) showed that mixing in the surface convec- 
tive zone of low-mass stars in the RGB phase is extended below 
the classical convective limit and that it is very fast compared with 
the nuclear time scales of either the hydrogen-burning shell or the 
^He-burning shell. Their model estimated that as the H-shell burns 
outwards the ^He will be destroyed. Although low-mass stars do 
produce considerable amounts of ^He on the main sequence phase, 
this will be all destroyed by the substantially deeper mixing they 
expect to occur in the red giant branch phase. These same results 
estimate that while 90% of the ^He is destroyed in IM© stars, only 
40-60% is destroyed in a 2Mo star model, depending on the speed 
of mixing. For stars above 4Mo, the bottom of the convective enve- 
lope is hot enough for nuclear burning (hot bottom burning). This 
leaves a mass range around 2Mo which is still expected to con- 
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tribute to ^He enrichment. Recent observations of Smiljanic et al. 
P009| ) have further compUcated the issue: they found a few 2.5- 
4Mo stars with low ^^C/^^C ratios, where CN-processed material 
has been brought from the H- shell to the envelope. Extra mixing 
mechanisms do not appear to operate in these stars, and additional 
mechanisms are assumed. Whether these would affect ^He is not 
known. 

The abundance of ^He can only be derived from measurements 
of the hyperfine transition of ^He^ which has a rest frequency of 
8.665 GHz. Detecting ^He^ in PNe challenges the sensitivity limits 
of all existing radio telescopes. The detection of ^He in any PN will 
require ^He/H~10~^, which is the abundance predicted by standard 
stellar models. After an intense search of ^He in PNe, the results 
are still unclear. [Bania et alTi ( t2010.) observed a sample of 12 PNe. 
^He was detected in only 2 of them. NGC 3242 was observed with 
Effelsberg a 100m dish from the Max Planck Institute for Radio 
Astronomy (MPIfR) and the NRAO 140-foot telescopes, however 
the observations are inconsistent with each other, implying a ^He/H 
abundance that is ~25% of the 10""^ abundance predicted. For the 
case of J320, ^He was detected at a 4cr level with the NRAO VLA. 
Composite ^He+ average spectrum for 6 PNe (NGC 3242, NGC 
6543, NGC 6720, NGC 7009, NGC 7662, and IC 289), using Ef- 
felsberg, Arecibo and GBT observations, consistently show ^He^ 
emission at the ~lmK level. 

^He is created in low mass stars and ejected during the late 
phases of their evolution. Galactic models predict a rapid increase 
over the past 5Gyr. This is however in contradiction to the observed 
abundances in the ISM. The cause is believed to be destruction of 
^He during the RGB phase, by hot bottom burning at higher masses 
and thermohaline mixing at low masses, leaving only a narrow 
mass range between. In this paper we present upper limits estimates 
for 3 PNe. The description of the observations made using the VLA 
in Section 2, the abundance calculations and results are presented 
in Section 3, the stellar models are presented in the discussion in 
Section 4 and the main conclusions are in Section 5. 



2 OBSERVATIONS USING THE VLA 

Motivated by this problem, three planetary nebulae were observed 
using the national radioastronomy observatory (NRAO) very large 
array {VLA) in D configuration. In Table [T] the targets and the pa- 
rameters of the observations are presented. 

The data was edited and reduced using the Common Astron- 
omy Software Applications (CAS A) package. The calibration tech- 
nique was slightly diff'erent for every object due to diff'erent S/N 
ratio in the phase and bandpass calibrators. 

• IC 418. For this source the phase calibrator (B0609-157) was 
used as a bandpass calibrator. The integration time was averaged 
every hour to get the bandpass solutions. The calibrator BOl 37+33 1 
was used to calibrate the intensity scale. The continuum intensity 
of this source was sufficient to perform the standard self-calibration 
technique. The self-calibration algorithm was applied to the contin- 
uum after the standard phase and bandpass calibration. 

Three loops of self-calibration were applied, the first 2 were only 
in phase, using the complete time interval for the first one and then a 
solutions interval of 300s. After this a self-calibration in amplitude 
and phase was applied using 150s for the solutions interval. Once 
this was obtained these solutions were applied to the line data, the 
continuum was subtracted and the clean image was created. 

• NGC 6572. For this source the phase calibrator (B 175 1+096) 



was used as a bandpass calibrator; the integration time was av- 
eraged every hour to get the bandpass solutions. The calibrator 
B 133 1+305 was used to calibrate the intensity scale. The contin- 
uum intensity of this source was sufficient to perform the standard 
self-calibration technique. The self-calibration algorithm was ap- 
plied after the standard phase and bandpass calibration. Several 
loops of self-calibration were applied. The first 4 were only in 
phase, changing the solutions interval from using the complete in- 
tegration time up to intervals of 21s. After this a self-calibration 
in amplitude and phase was applied three times, also changing the 
intervals of the solutions, starting using 150s. Ultimately the best 
result for the self-calibration was obtained using a solution interval 
of 21s in phase only and a solution interval of 30s for the amplitude 
and phase. Once this was obtained these solutions were applied to 
the line data, the continuum was subtracted and the clean image 
was created. 

• NGC 7009. For this source also the phase calibrator 
(B21 36+006) was used as a bandpass calibrator; the integration 
time was averaged every hour to get the bandpass solutions. The 
same calibrator was used to calibrate the intensity scale. For this 
data set, a flux calibrator was not observed. But as mentioned in 
Table [T] there were more observations taken for these sources. In 
the case of NGC 7009 a day before of the observations analysed 
here a flux calibrator was observed; this was used to calibrate the 
phase calibrator (same for both observations) and the flux obtained 
was imposed into the phase calibrator for the data set used here. 
Therefore the same calibrator was used as a phase, bandpass and 
amplitude calibrator. For this source a self-calibration technique 
was applied, but it did not improve the data. Therefore the image 
and spectrum presented here were made using only the standard 
calibration. Once this was obtained these solutions were applied to 
the line data, the continuum was subtracted and the clean image 
was created. 

The continuum images are completely limited by the finite 
dynamic range of the strong continuum emission and thus are not 
limited by thermal fluctuations. The line data, after subtraction of 
the continuum, are limited by thermal noise. This explains why the 
continuum rms values are larger than the line data values (see Table 



3 RESULTS 

The result of the observations are presented in Figure [T] showing 
the emission cube with the continuum subtracted. As a reference 
the emission of the continuum is over plotted in red contours. The 
spectrum presented on the right of every image is the emission com- 
ing from the centre of the contours, where the emission of the line 
should be. The ^He^ line was not detected for any of the objects, 
therefore only upper limits for the density values were calculated. 
We used a 3cr level of the rms as an upper limit for the three sources. 
The ^He^ column density can be obtained using: 



7V(^He+) : 



gi + gu ^nkv^ 



hc^A, 



' r 

- Tb(v)Jv, 

ul J 



(1) 



where gu =1 is the g value or magnetic moment of the upper 
state, gi =3 is the g valu e or magnetic moment of the lower state, 
Aui =1.95436xl0-i2s-i ( |Gould|l994 t, and Tb{v) is the brightness 
temperature profile of the line. For a gaussian line profile. 



J TB(v)dv=\/2,^T'sAv, 



(2) 



^He: Does the problem persists? 3 




-6x1 



30^ 29^ 28' 27' 25' 
J2000 Right Ascension 





1.5x10"^ 



08= 07= 06= 05= 04= 
J2000 Right Ascension 




20 40 
V_LSR (km/s) 




14^ 12= 10= 08= 06= 04= 
J2000 Right Ascension 




-40 -20 

V_LSR (km/s) 



Figure 1. VIA observations of IC 418 (top), NGC 6572 (middle), and NGC 7009 (bottom). On the left side the image shows the line emission of the source, 
continuum subtracted. The red contours in the centre show the position of the continuum emission for each source. On the right side is the spectrum at the 
peak continuum position of the source. The ^He+ line should be at 62 km s~^ for IC 418, at -8.7 km s~^ for NGC 6572 and at -46.6 km s~^ for NGC 7009. 
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Table 1. VLA Observations Parameters. 



Parameters IC418 NGC 6572 NGC 7009 



Date* 


15-12-2009 


13-12-2009 


26-12-2009 


Total observed time (hr) 


6 


6.2 


4.6 


Configuration 


D 


D 


D 


R.A. (J2000) 


05:27:28.20 


18:12:06.36 


21:04:10.8 


DEC (J2000) 


-12:41:50.26 


-^06:51:13.01 


-11:21:48.25 


FWHM of primary beam (arcmin) 


5.4 


5.4 


5.4 


Size of the PNe (arcsecs) 


30x6 


18.9x6.9 


12x7 


FWHM of synthesised beam (arcsec) 


13 


9 


9 


Total bandwidth (MHz) 


6.152 


6.152 


6.152 


Number of channels 


62 


62 


62 


LSR central velocity (kms~^) 


62.0 


-8.7 


-46.6 


Observed v (MHz) 


8662.6338 


8665.6132 


8666.0602 


Spectral resolution (kms~^) 


3.378 


3.378 


3.378 


Flux density calibrator 


0137-^331 


1331-^305 


1733-130 


Phase calibrator 


0609-157 


1751-^096 


2136-^006 


Line channel rms (mJy beam~^)** 


1 


0.4 


0.8 


Continuum rms (mJy beam~^ ) 


2.0 


1.53 


1.3 



*More observations were obtained, due to the lack of calibrators or bad data they were not considered here. 
Continuum subtracted 



where is the brightness temperature at the centre of the line and 
Av is the full width at half power. The brightness temperature Tg is 
related to the observed beam-averaged brightness temperature Ti, 
by ( |Rood etal.|1995| l 



Tl- 



01 



(3) 



where 6}, and 6s are the beam and source angular radii. 

The number of ^He^ atoms per unit volume, ;i(^He^), can be 
obtained by dividing the column density A/^(^He^) by the averaged 
optical path, < A5 >, through the source. Representing a FN as a 
homogeneous sphere of radius R^^^^D, where D is the distance, the 
optical path at a position angle 6 is I^s{6) - 2 ^]B} - (OD)^, and 
the optical path averaged on the source results < A^ >= ttR/I = 
nOsD/l. 

The expression for ;i(^He^) is then 



w(^He+) : 



n g i + g, kv' 61 + 6; 

gu hc^Aui ml 



TlAv. 



(4) 



Rearranging Equation |4] and using the numerical values, the 
final expression can be presented as follows. 



w(^He+) 



22.7 



Av 

\mK/ \kms~^ 

3 / /) x2 



D 

kpc 



(5) 



\ arcsec/ \AV' I \ el) 



Table [2] shows the calculated ^He^ density using Equation [s] 
for the three sources. The distance values used were L2kpc for 
NGC 7009 ( |Galli et al.|19"97]) 1.2kpc for NGC 65 72 pajian et al.| 
[1995^ and l.Skpc for IC 418 ^Guzman et al.|2009l ). To calculate the 
abundance we need to divide by the density. The values pre- 
sented in Table [2] were taken from Galli et al.|( |1997| ); |Balser et al.| 
([2006};[Morisset & Ge orgiev ( 2009 ) for NGC 7009, NGC 6572 
and IC 418 respectively. The abundance is derived from the 
thermal free-free bremsstrahlung radio continuum emission, and is 
proportional to the square of the density. To obtain these values, 
the PNe where modelled using a numerical code, called NEBULA 



for the case of NGC 6572 and NGC 7009. This code calculates the 
radio recombination line and continuum emission from an ionised 
plasma. The modelled nebulae are composed solely of hydrogen 
and helium. The numerical grid uses a three-dimensional Cartesian 
coordinate system that allows for the specification of a completely 
arbitrary density and ionisation structure. For the case of IC 418 its 
abundance was modelled using Cloudy _3D ( Mori sset|2006| ). 
Once we calculated the ^He^/H^ abundance, we can compare 
the values with the theoretical calculations and the stellar models 
predictions. 

Figure [2] shows the ^He abundance of the three PNe from this 
work, they are presented using the red symbols. The mass values 
used are 1.7+0.3 for IC 418 (Moriss et & Georgiev|2009j ), 1.0+0.2 
for NGC 6572 (|Hyung et al.„1994j ) and 1.4+0.2 f or NGC 7009 

1 Galli et al.||1997| ). For comparison the PNe from [Balser et aL| 
1997 ) are also presented (using green symbols), where the mass 
estimates are from G alli et al.|(l997 ). All these measurements rep- 
resent upper limits points. For all the ^He^/H+ values calculated 
and taken from the literatu re it has been assumed that ^He/H^ 
^HeVH^ ( |Balseretal.|l994t . 

PNe IC 418 and NGC 7009 have values consistent with Pop- 
ulation I and II theoretical predictions. For IC 418, this is the first 
upper limit calculation and it agrees with the stellar model predic- 
tions. For NGC 7009 a ^He/H abundance of < 8.24 x 10"^ was 
calculated. Previous observations of this source gave estimates of 
< (5.8 - 31) x 10~^, therefore our estimated value does not improve 
earlier estimates. In the case of NGC 6572 the abundance calcu- 
lated is lower than the theoretical models predict. An upper limit 
for this PNe was calculated by Bals er et al.| ( [2006] ): the value they 
estimate is < 10~^. We found an upper limit of < 2.4 x 10~^. 



4 DISCUSSION 

The stellar evolution models are also presented in Figure [2] The 
curves labelled as Pop I and II show the standard abundance of ^He 
taken from Weiss et al.|(T996j ). The Big Bang line is the primordial 
value of ^He/H. 

In order to understand the chemical evolution of the Galaxy, it 
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Table 2. Parameters used to calculate the ^He abundance in the PNe. 



PNe 


D 


M 


Tl 


Av 


^(H+) 


«(3He+) 


3He+/H+ 




(kpc) 


(Mo) 


(mK) 


(kms-i) 


(cm-3) 


(cm-3) 


(10-3) 


IC418 


1.3 


1.7+0.3^ 


27 


30 


3x10^ 


<12.99 


<4.33 


NGC 6572 


1.2 


1.0±0.22 


150 


15 


22x1 0^ 


<5.38 


<0.24 


NGC 7009 


1.2 


1.4+0.2^ 


465 


43 


3.62x10^ 


<29.85 


<8.24 



^Morisset & Georgiev|j2009) , ^Hyung et al.]il994t , ^GaUi et al.|fl997t 



is crucial we understand the contribution of PNe to the ^He abun- 
dance. Chemical evolution models started including only processed 
primordial deuterium (D) as a source of ^He synthesis in low mass 
stars. But even without including any stellar producti on of ^He , 
these models showed that ^He/H increases with time (Vangioni 



red giant stars that are descendants of Ap stars does reconcile the 
measurements of ^He/H in Galactic H ii regions with high values of 
3 He in a few PNe. 



Flam et aLl|l994 l. Because the abundance of ^He in Hii regions, 



the ISM and the proto- solar system is much lower than predicted, 
new chemical models were proposed. 

These models include non convective mixing processes that 
either prevent the buildup of ^He on the main- sequence or destroys 
3 He along the upper red giant branch ( [Boothroyd & Sackmann 
[l999||Lagarde et al.|2012| ). 

The blue curve labelled ExtraMix represents the results of stel- 
lar nucleosynthesis calculations using deep mixing by [Boothroyd &| 
|Sackmann| ( |1999| ). The detailed calculations indicate that the mass 
dependence of the destruction of ^He is very strong, and the result- 
ing abundance decreases sharply for masses below ~2.5Mo. It is 
clearly seen in this Figure that the observed PNe have not suffered 
any depletion. Almost all of the PNe are more than a factor of ~2 
above the non-standard curve (the Extramix predictions). 

One suggestion to solve this is by adding some additional mix- 
ing mechanism during the first ascent of the RGB. By allowing 
material to circulate below the convective envelope after the first 
dredge-up has occurred, the surface abundances can be modified as 
the star approaches the tip of the RGB. Observationally, the exis- 
tence of this no n-convective mixing is well established (e.g. |Grat^ 
[ton et al.|2600t as carbon, lithium and the ^^C/^^C ratio are all ob- 
served to change as stars reach the tip of the RGB. That this pro- 
cess should also aff'ect the ^He abundance seems clear as to aff'ect 
these elements and isotopes one must reach temperatures where 
^He burning reactions will occur. However, the physical cause of 
this non-convective mixing has proven elusive and many possible 
mechanisms, such as rotation (e.g. |Palacios et al.|2006| ) have been 
suggested. 

The hydrodynamic simulations of [Eggleton et al.| ( |2006| ) 
pointed to a particularly convenient process for this extra mixing. 
Their simulations showed that a mean molecular weight inversion 
caused by the destruction of ^He could drive circulation of material. 
This process (subsequently identified with thermohaline convection 
by |Charbonnel & Zahn|2007^ would allow low-mass stars to de- 
stroy the 3 He that they produce during the main sequence phase. 

Newer stellar models from| Lagarde et al.] (|2012l) analyse the 
evolution of ^He and its destruction due to the eff'ects of thermo- 
haline instability and rotation-induced mixing. Their models show 
that thermohaline mixing can induce significant depletion of ^He in 
low- and intermediate-mass stars. Although these stars remain the 
net producers of ^He, their contribution to the Galactic evolution 
of this element is highly reduced compared to classical theory. It is 
claimed in their paper that thermohaline mixing is the only physi- 
cal mechanism known so far able to solve the so-called "The ^He 
Problem". Importantly, its inhibition by a fossil magnetic field in 



We have used the STARS stellar evolution code (Eggl etonj 
[T97T] [Pols et al.||T995l [Stancliff^e & Eldridge ,2009 ) to produce a 
set of stellar models in the mass range 1-3.5 solar masses at inter- 
vals of 0.25 solar masses and with a metallicity of Z=0.02. These 
models have been computed from the pre-main sequence to the first 
thermal pulse on the asymptotic giant branch. Thermohaline mix- 
ing has been included at all stages of the evolution using the dif- 
fusive formahsm of |Ulrich|(T972l) and [Kip penhah n et al.| ( [T980l ). 
Following the work of Ch arbonnel & Zahn| ( |2007a| , we set the free 
parameter of the formalism to Q = 1000, a value which not only 
reproduces the compositional changes in field giants, but which has 
also been shown to reproduce the diff'erence in behaviour between 
C-normal and C-rich metal-poor stars ( [Stancliff'e et al.| 



2009t 



well as the abundance changes in globular cluster stars ( Angelou| 
|et al.|201l] ). Note that this value is not supported by current hydro- 
dynamical simulations of thermohaline mixing ( |Traxler etal.|201l( 
[Denis senkov & Merry field|201 If and that some other process (or 
processes) could be the cause of abundance changes on the giant 
branch (e.g. magnetic fields, jCharbon nel & Zahn'2007bVNord haus[ 
[et al.|f2 008 ). These stellar models using thermohaline mixing are 
also represented in Figure [2| (light-blue line). 

As shown in Figure [2]the ^He abundance of two PNe (NGC 
6572 and NGC 6720) fall below the line that the model predicts, 
therefore at least for these two PNe the thermohaline mixing model 
fails to reproduce their ^He abundances. 



5 CONCLUSIONS 

Upper- limits for the ^He abundance of three PNe have been esti- 
mated using the ^He^ emission line observed with the VLA towards 
the PNe IC 418, NGC 6572 and NGC 7009. In total one strong de- 
tection (PN J320) and seven upper limits of the ^He abundance on 
PNe have been reported in the literature, including this work, and 
are consistent with standard stellar evolution calculations. These 
estimates are two orders of magnitude larger than what is found in 
Hii regions ( [Balser et al.|1997j l, the local ISM ( jGloeckler & Geiss[ 
[1996[ ) and the proto-solar system { [Geiss|1993] ). 

The resulting evolution of ^He is only consistent with the val- 
ues determined in pre- solar material and the ISM if more than 96% 
of the whole population of stars with mass below ~2.5Mo have un- 
dergone enhanced ^He depletion. These results suggest that more 
observations are needed in order to make a strong conclusion about 
the stellar evolution models. 
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PN6(Balseretal.,1997) ngctooq 

: ~r 




1 1.5 2 

stellar Mass (Msun) 

Figure 2. Abundance of -^He/H versus main-sequence masses. The red sym- 
bols represent the work from this paper. The horizontal bar represents the 
mass range of the progenitor star. The green symbols represent the sam- 
ple of 6 PNe (N GC 3242, NGC 6543 , NGC 6720, NGC 7009, NGC 7662 
and IC 289) from jBalser et al.|fl997) , using the mass estimates from |GallI] 
|et al.|jl99 7L The curves as Pop I and II show the standard abundance of 
^He taken from Weiss et al. (1996 1. The blue curve labelled ExtraMix rep- 
resents the results of stellar nucleosynthesis calculations using deep mixing 
by ,Boothroyd & Sackmann { 1999) . The light-blue curve represent the stel- 
lar models using thermohaline mixing. The Big Bang line is the primordial 
value of ^He/H. 
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